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1 Introduction 



A basic tool of quantum field theory is the generating functional for multi-time correla- 
tion functions. In this paper, we show how the time-dependent variational principle of 
Balian and Veneroni |IJ can be implemented in the context of quantum field theory out 
of equilibrium. From this approach, we derive approximate non perturbative dynamical 
equations for expectation values and correlation functions of local and composite opera- 
tors. These evolution equations could be useful to evaluate the fluctuations of the matter 
density in cosmological scenarios @. They could also provide informations about the 
non-equilibrium dynamics during the expansion of the quark-gluon plasma 0, [3J. We 
consider here the case of a self-interacting scalar field in a Minkowski metric described by 
the hamiltonian H = J d d x7i(x) with 

1 1 -» Tn^ h 

H(x) = - n 2 (f) + - (v$(f)) 2 + -f $ 2 (f) + - $ 4 (f) . (l.i) 

We work in d spatial dimensions. The constants ttlq and b are respectively the bare mass 
and the bare coupling constant. 

2 Variational evaluation of the generating functional 

Let us define an operator A(J, t), functional of the sources J and function of the time t, 
having the form : 

A(J,t) =TexpiJ™dt' y d d x J*(x,t')<$> H (x,t',t) + J n (x,t')U H (x,t',t) 

+ J d d x d d y J* n (f, y, t') (§ H (x, f , t)U H (y, t', t) + U H (y, t', t)$ H (x, t', tj) , 
+J**(x,y,t') <S> H (x,t',t)$ H (y,t',t) + J uu (x,y,t>) U H (x,t',t)U H (y,t',t)} 

(2.1) 

where T is the time-ordering operator and $ H (x, t', t) and U H (x, t', t) are respectively the 
field and momentum operators defined in the Heisenberg representation. They satisfy the 
boundary conditions : $> H (x,t,t) = $(x) and U H (x,t,t) = H(x), $(x) and Tl(x) being 
the operators in the Schrodinger representation. The local as well as bilocal sources J{t') 
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are turned on between t' — t and t' = +00. The generating functional for the causal 
Green functions writes 

Z(J,t )=Tr(D(to)A(J,to)) , (2.2) 

D(to) being the initial state. The statistical operator D(t ) may represent a thermal 
equilibrium or a non-equilibrium situation. At zero temperature, it reduces to a projection 
operator. Besides the sources J, the generating functional Z depends on the initial time 
t . We want to evaluate the functional derivatives of 

W(J,t ) = -i lnZ(J,t ) , (2.3) 

since W( J, t ) is the generating functional for the connected Green functions. Its expan- 
sion in powers of the sources writes : 

W(J, t ) = -in + dt' { J d d x J$(x, t') ip(x, t') + J d d x d d y </$$(£, y, t') G(x, y, t') + ...} 

+ % - J J^°° dt' dt" { J d d x x d d x 2 J*(fi, ,U(x 2 , t") C| 4 (fi, x 2 , t', t") + ...} 

+ ... 

(2.4) 

where n = TrD(t ) ( we choose not to normalize D(t )). We use the following notations 
for the expectation values of one and two field operators : 

( p(Z,t) = —Tr($ H (x,t,t )D(toj) , (2.5) 
no x ' 

G(x, y,t) = — Tr ($ H {x, t, t ) $ H (y, t, t ) £>(*„)) - <p{x, t) <p(&, t) . (2.6) 
no 

The two-time causal functions between two field operators are defined according to : 

C| $ (fi, x 2 , t u t 2 ) = —Tr(T ^ H (x u t u to) <5> H (x 2 , t 2 , t ) D(t )) - <p(x u h) <p(x 2 , t 2 ) . 

(2.7) 

Similarly, we define C| n , C^ u . We define also the three-point and four-point two-time 
correlation function C 3 and C 4 according to : 

^ 5Mx, V)5J$*(y, z, W) 1 J=0 = C3(f ' ^' r) 

+V(y,t") (C\x,z,t',t») + V(x,t>) ^z,t»)) + V&t") (C 2 (x,y,t',t") + ^,t')^y,t")) , 

(2.8) 



8^W 

5J$$(x, y, t')5J^(z, u, i 



__ = C 4 (x,y,z,u, t',f 



+ (C\x , u, e, t" ) + , t" )) (C 2 (y, z, if, f ) + ^(y, p(if, f ' ) 



(2.9) 



We will use the time-dependent variational principle of Balian and Veneroni |l] to ob- 
tain an approximation for our quantity of interest, Tr(D(to)A(J, to)). We define therefore 
the action functional 

Z(A(t),V(t)) = Tr(A(t )V(t )) + Z dyn , (2.10) 

with 

Z>dyn = Tr 



J t °°dtV(t) A(t),J d d xH(x) +iA{t) [YjJ^qM , (2.11) 

where we have written in a compact form the term which involves the sources J : 

Jj(t)Qj = J d d x £)$(£) + J n (x, t)Tl(x)) 

+ J d d x d d y J* n (f , y, t) ($(f)n(y) + n(y)$(f)) (2.12) 
+ J d d xd d y (J* $ (f,y,t)$(f)$(y) + J nu (x,y,t)Tl(x)Tl(y)) ■ 

The variational quantities of the functional Z, which depends on the sources J, are the 
observable-like and density-like operators A(t) and T>(t). We look for the stationarity 
of Z with respect to variations of A(t) and T>(t) subject to the boundary conditions : 
A(t = +oo) = 1 and T>(to) = D(to), where D(to) is the initial statistical operator which 
we assume to be given and equal to a Gaussian density matrix. The generating functional 
for the connected Green functions will be approximated by W(J,t ) = —i lnZ st . 

We restrict ourselves to trial operators A(t) and T>(t) which are exponentials of 
quadratic and linear forms of the field operators $(x) and U(x) (which we shall loosely call 
Gaussian operators). Let us introduce the operators % = A(t) T>(t) and T c = V{t) A{t) 
||. We define the following expectation values : 

n(t) = TrT b = TrT c , (2.13) 

Vb {x,t) = -^Tr{T b §{x)) , (2.14) 
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n b (x,t) = ^-Tr(T b U(x)) , (2.15) 
n(t) 

G 6 (£,j?,*) = ^ Tr{T b ${x)®{$)) , (2.16) 

S b (x, y,t) = ^- Tr(T b IL(x)IL(y)) , (2.17) 
n(t) 

T b (x, V,t) = 4x Tr(%($(x)fi($) + fi{y)®{x ))) , (2.18) 

with $(x) = $(x) — ip b {x,t) and n(x) = II (a?) — 7r b (x,t). Similarly we define the ex- 
pectation values associated with T>(t)A(t) and indexed by c by replacing T b by T c in 
( |2.14D -( ^l8l ). For £>(£) = 1, one has p b (x,t) = (p c (x,t) = Lp a (x,t) (and similar defini- 
tions for 7r a , G a , S a , T a ). For A{t) = 1, one has f b (x,t) = (p c (x,t) = ipd(x,t) (and similar 
definitions for % d , Sd, T d ). A Gaussian operator is completely characterized by a set 
{tp, 7r, G, S, T}. It will be convenient to introduce a vector a(x, t) 

(p(x, t) 

a(x,t) = I I , (2.19) 



-Z 7r(x, t) 



and a matrix S(5;, t) 



2G(x,y,t) -iT(x,y,t) 
-iT(y,x,t) -2S(x,y,t) 
We give in appendix the expressions of c^, a c , S 6 and S c as functions of the independent 
variational quantities a a ,S a and otdj^d which characterize A and T> respectively. 

With this choice of the trial spaces for A(t) and T>(t) , the Wick theorem allows us to 
express the functional Z in the form : 



Z = n(t ) + / dt 



-OO 



to 



dfi f 

-jj;\v(t)=cte - i n(t) J d d x (£ c (x , t) - £ b (x, t)) 



(2.21) 



+in(t) J d d x\d d X2 JC c (xi, x 2 , t) 
where the energy density S(x, t) = Tr (T>(t) 7i(x, t)) is given, for the Hamiltonian density 
(OD, by : 

£{x, t) = l - n 2 (x, t) + \ (Vif) 2 + ^ cp 2 (x , i) + A , t ) 

+1 S {x, Z,t)-~ AG(x, y, t)\ x=y + ^ G(x , f , t) (2.22) 
+- G 2 (f , f , t) + - (p 2 (x, t) G(x, x, t) . 
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The last term of the functional Q2.21]) , K c = Tr (T c J2j Jj(t) Qj^j = J d d Xid d x 2 }C c (xi, x 2 , t), 
involves the sources J and the expectation values indexed by c. 



3 Dynamical equations 



By varying the expression (|2.21|) with respect to n d , a d , E d , with the boundary conditions 



n d (t ) = n , a d (t ) = a , E d (t ) = S , (3.1) 

where no,«o et So characterize the initial Gaussian state D(to), we obtain the evolution 
equations for n a ,a a and E a . Integrating ( |2.21| ) by parts and varying with respect to 



n a ,a a ,E a with the boundary conditions 

n a (t = +oo) = 1 , a a (t = +oo) = , ~-\t = +oo) = , (3.2) 

we obtain the evolution equations for n d , a d and E d . In general the evolution equations 
for n d , a d , E d and those for n a , a a , S a are coupled. The solutions n d , a d , E d and n a , a a , E a 
depend on the sources. 

By combining the evolution equations for n d ,a d , E d with those for n a , a a , E a , the dy- 
namical equations for the expectation values at,, and a c , S c can then be written in the 
following compact form : 

i 6ib = t w b , (3.3) 
i a c = T (w c - w c K ) , (3.4) 
iZ b = 2(z b H b T-TH b Z b ) , (3.5) 
i ~ c = 2 (~ c (H c - T K ) t-t(H c - T k ) S c ) , (3.6) 
where r is the 2x2 matrix 



(3.7) 



1 

-1 

The vector w b and the matrix T~C b are defined from the variation of < H > b = Tr(T b H) : 



5 < H > b = J d d x w b (x, t) 5p b (x, t) - - J d d x d d y trH b (x, y, t) SE b (x, y, t) . (3.8) 



Similarly, w c ,7i c and -u^Z^- are defined from the variation of < H > c = Tr(T c H) and 
from the variation of the source term K c . 

Therefore, using a convenient parametrization for the two variational objects T>(t) 
and A(t), we have been able to obtain dynamical equations in a very compact form even 
for finite values of the sources. The dynamical equations for a>d, and a a , S a have a 
more complicated form. They are given in appendix B, where we write also the explicit 
expressions of the vector w and the matrices 7i and X C K in the case of the A$ 4 theory. 
In spite of their form, the equations ( p.3|) -( pl6|) are coupled because the solutions a?,, E& 
and a c , S c do not satisfy simple boundary conditions. The expansion in powers of the 
sources of the stationarity conditions ( |3.3|) and (|3.6|) will provide approximate dynamical 
equations for the expectations and correlations functions defined in eqs. ( p.5|) -( ^9|) . 



4 Expansion in powers of the sources 

We will use upper index (°) for the solutions of the dynamical equations when there are 
no sources. The limit with vanishing source corresponds to = and S^ l(0> = 0, 
= a (0) = 4°) a nd s£ 0) = = sf; we have also w b ^ = w c ^ = w<® and 
n b(o) = n do) = ^(0) In this limit; the dynamical equations (^TT|) and ( |B^) for a d et E d 
become (we supress the index d) : 

m'^r/, (4.1) 

i S<°) = - [(S® + r) ft(°> (E® - r) - (»(°) - r) H<® (z {0) + r)} . (4.2) 

These equations are the analog for the A<3> 4 theory of the time-dependent Hartree-Bogoliubov 
(TDHB) equations in non-relativistic physics. They are equivalent to the dynamical equa- 
tions obtained in reference || where the authors used an alternative form of the Balian- 
Veneroni variational principle suited to the evaluation of single-time expectation values 

The first derivatives of W(J, t ) = —i \nZ st with respect to the sources are equal to 
the expectations values with the index c. Indeed, the functional Z depend on the sources 



both explicitly and implicitly since the approximate solutions for V{t) and A(t) depend 
on the sources. However at the stationary point, only the explicit dependence contributes 
to the first derivative : 

-^L = iTr(V(t)A(t)Q j ) , (4.3) 

which gives for instance : SJ 5 J^ ^ = ip c (x,t). The expressions for the second derivatives 
of W are much more complicated. The introduction of sources coupled to the composite 
operators $(?/), &(x)U(y) and U(x)U(y) together with eq. (|4.3| ), allows us to obtain 
dynamical equations for two-time correlation functions with three or four field operators 
merely from the expansion of the expectation values a c et H c at the first order in powers 
of the sources. 

From the first order corrections a c — and S c — S^, we define the two-time corre- 
lation functions (3 and £ (j3 is a vector and £ is a matrix) : 

a c (x, t)- a (0) (x, t)~i l°° dt" y d d y /?*(£, y, t, t") J (y, t") + /9 n J n 

+ J d d Xl d d x 2 (3 m (x, xi, x 2 , t, t") J**(f x , f 2 , t") + /?* n J* n + /3 nn J nn | 

(4.4) 

/•+oo ( /■ 

E c (x,y,t)- E(°\x,y,t)~i jf dt" £*(£ , y, , f l5 i, t") J $ (x l5 t") + S n J n 

+ 1 d d Xl d d x 2 , jf, fi, f 2 , t, t") J $<I> (£i, f 2 , i") + S $n J* n + E nn J nn | 

(4.5) 

We have 

(3?*(x,y,z,t,n = Z* 1 (x,y,z,t,t») . (4.6) 

The functions (3 and £ provide approximations for the exact two-time correlation functions 
C 2 , C 3 and C 4 defined by eqs. fl2?7])-Q: 



C£*(£,jf,t,0^#(z,£M'') , ( 4 -7) 

C 3 (f,y,f,M" )~ 0**(x,frz,t,V) 

-ipWfaf) (pf(x,z,t,r) + ^°\x,t)^ \zX)) , (4-8) 
-^°) (F , f ' ) (/5f (x , £ t, t" ) + <pV> (x, t)<p<® (y, t" 
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C A (x , y, z,u,t, t" ) ~ Eff (f , y, z, u, t, t" ) 

- (bf (x , z, t, + , t) ^ (0) (^ «**)) (6f (3/, 3, t, + V (0) (y, *) V (0) (^ 

(4.9) 



The expansion in powers of the sources of the dynamical equations ( 3.4j) and (|3.6| ) for 
j3 c et E c yields : 

iSj3 c = r (5w c - 5w c K ) , (4.10) 

L . (4.11) 

+ ~ {0) (^ c - 51 C K ) t-t (5H C - 5X C K ) H (0) 

In these equations the matrix ri has to be evaluated for the TDHB solutions a^°\S^ 

ofeqs. (PD-flPQ. The variations of w c K (x,t) (eqs. ( P^ )- (P35|) ) and T c K (x,y,t) (eqs. 

( B.16|) -( [B.18| )) with respect to </(£") will give terms proportional to 5(t — £"). From the 

variations 5w c and 6TC C , we define the matrices t^, T^, r ij,k, Rij,ki (which are the analogs 

of the RPA kernel of Balian and Veneroni Hill) : 



1 



8w-(x, t) = tij(x, t) 5a c j(x, t)-- d y T iJk (x, y, t) SE c kj (y, x, t) , 



(4.12) 



1 



5H c ij {x, y, t) = r ijtk (x, y, t) 5a c k {y, t) - ~ \ d d z R ijM {x, z, t) 5E c lk (z, y, t) . 



(4.13) 

These matrices depend on the TDHB solutions and H 1 - ). Their expressions in the 
case of the A$ 4 theory are given in appendix B. With these notations, the dynamical 
equations for the two-time and three-point correlation functions /?* and E* write : 



i£p*(x,y,t,t")=r 



t(x, t) f3 (x, y, t, t") - - T (jk (x, x, t) Ejyj(x, x, y, t, t") 



+i 



6{x-y)5{t-t") 



E* H (0) t-t E* 



+2 (0) r (fc 0* t - t r (fc (3% S<°> 



, (4.14) 



K ) 1 ~ i • ,(* ^it) ' • (4-15) 

z z J 2j 

We obtain similar equations for the three-point functions /? <1> *,/3 <1>n , /3 nn and the four- 
point functions £ <I> *,E <I>ri e t E nn . These dynamical equations are not sufficient since 



the boundary conditions are a a (t = +00) = 0,S a ; 1 (t = +00) = and ad(x,t ) = 
cxq(x), Ed(x,y,to) = Eo(x,y). We will use the dynamical equations for the functions / 
defined from the expansion of S" 1 at the first order : 

E- 1 (x, y,t)~ i J t °° dt" I J d d Xl f (£, y, f 1 , t" , t) J* (f ! ,t") + l n J n 

+ J d d Xl d d x 2 y, x u x 2 , t" , t) J**(2 U 2 2 , t") + Z* n J* n + Z nn J nn | 

(4.16) 

We have also at the first order : 



a — a 



and 



At the initial time to, 



and 



a c (t ) - «o - (t S a 1 - H a 1 S 



(4.17) 
(4.18) 
(4.19) 
(4.20) 



E c (t ) - E - 

This yields the following relations between the functions /3*(t ,£"), S*(t ,t") and the 
function /*(£" , t ) : 

£i,*o,*") = / d d z 1 d d z 2 (r f(x, z 1: xi, t" ,t ) <K-?l - ^2) - Z*(f, ^i,xi,t",t ) 2 (Fi,z 2 )) "0(^2) 

(4.21) 



S*(f,y,fi,t ,t") = - j d z x d z 2 (E (x,zi) ~ t) 

x (So(z!,^) l'^(z 2 ,y,x 1 ,t , \t ) + 5(z 2 - y) l*(z u y, 2 X , t" , t ) r) 

(4.22) 

We have similar relations between the functions and Z**. 

To lowest order, the dynamical equation for S" 1 reads : 



d 



H c -H b -T K -2 E^tH^ + 2 H^tE- 1 



(4.23) 



From the expansion ( 4.16 ) for S a x and the variations £7Y C — £7Y 6 , we obtain 

% ±1(2, y, x x ,r , t)ij = -2 r y , fc r fci l lm 0$ - R i3M (~ (0) It-t E<® l) 

+ [-2 lrH {0) + 2 H {0) rl + i 5X C K S(t - t" 



Ik 



(4.24) 
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We have l(x, y, t" , t) = for t > t" and : 

l**{x,y,x 1 ,x 2 ,t",t") = l , l m {x,y,x u x 2 ,r,t") = 2i , l nu {x, y, x u x 2 , f , f ) = -1 , 

(4.25) 

/* (x , y, x x , t" ,n = l n (x, y , x x , i" , t" ) = . (4.26) 

Therefore, to obtain approximations for the two-time causal correlation functions 
C 2 ,C 3 and C 4 , we need first to solve the TDHB equations (fOD and fl4|) for a (0) and 
and to evaluate the matrices r,R,t and T. Then we solve the equation ( 4.24|) 



for I (x, y, x±, f , t) and I (x, y, Xi,X2,t" ,t) (t" > t' > t ) backward from t" to t 



with the boundary conditions (|4.25| )- (^4.26|) . From the relations (|4.21|) and ( [4.22 



we 



obtain @®(x, x±, t , t"), {3®®(x, X\, x 2 , t , f) , S*(a;, ^, Xi, f 0) ^") and E**(x, xi, x 2 , t , t") 
and with these boundary conditions we solve the dynamical equations for 0(t,t v ) and 
£(£, i") forward from t to t" . 

In the TDHB approximation, the three-point function C 3 and the four-point function 
C 4 vanish. Using the Balian-Veneroni variational approach, we are able to obtain ap- 
proximations for C 3 and C 4 given by eqs. (|4.8|) and ( [4.9D which differ from this naive 
result. 

In order to obtain approximations for the two-time anticausal functions, it is sufficient 
to replace in the functional ([2.10|) -( p.ll|) the term K c by the analog term Kj,. We can 



therefore derive a consistent approximate formula for the response functions which involve 
the retarded commutator. 

When the two times coincide in the function /3f , we will obtain an approximation for 
the function G(x, x±,t) which will differ from the TDHB solution (they satisfy however 
the same initial condition G(x, Xi,t ) = G^(x, Xi, t ) = G (x,xi) ). We investigate 
the static situation in a subsequent paper where we use a form of the Balian-Veneroni 
variational principle which optimizes both the initial state and the dynamics JIJ . We will 
examine the divergences which appear in the variational equations and study how the 
approximations obtained for the functions G, C 3 and C 4 are related to a few physical 
quantities of the theory, i. e. the renormalized mass of the particle and the renormalized 
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coupling constant. 

In conclusion, using an appropriate parametrization, we have shown how to evaluate 
variationally the generating functional for multi-time correlation functions in a quantum 
field theory out of equilibrium when the initial state is given and equal to a Gaussian 
state. We have obtained dynamical equations in a compact form whose expansion at the 
first order in power of the sources provides approximations for the two-time correlation 
functions. 
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Appendix A 

In this appendix we give the expressions of the expectation values a b , a c , E b and H c as 
functions of a d , a a , E d and S . 



a b = (E a - t) 




— ®d + (S d + r) - a a , 



(A.l) 



a c = - r 



) 




— a a + (S + r) - a d , 



(A.2) 



where r is the 2x2 matrix : 




(A.3) 



For the matrices S, we have the following relations : 



E b - r = (E a - t) - - (E d - t) 



(A.4) 



S c - t = (E d - r) - — — - (S a - r) . 



(A.5) 



Appendix B 
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In this appendix, we give the evolution equations for a d , E d and a a , E a . 



2i a d 



[E d + r) (H c - T K ) (E d - r) - (E d - r) H b (E d + r) 



r a 5 _ 



+E d (w c -w b )+r (w c + w b ) - (E d + r) u& 

S, + r) (H c -X^) - r) - (S d - r) H b (E d + r) 



2i a a = (E a - t) (H c - T K ) (E a + r) - (E a + r) H b (E a - r) 



+E a (w b -w c )+t (w c + w b ) - (E a - t) w~ k 

E a -r) (H c -T K )(E a + T)-(E a + T) H b (E a - r) 



The expressions for the vector w and the matrices H and X are the following 

w b (S, t), = S< g > b = -f\x, t) , w\x, t) 2 = i S< ?_ = t) , 

d(fi b (x,t) dir b (x,t) 

Ti"(x.y.0 7:i --2- 6 " 



5E b (x,y,t)ji ' 

H b (x, y, t) u = - = ^ y, , 

5G b (y,x,t) 2 

6T b {y,x,t) 

For a self-interacting scalar field, we have : 

f b (x, t) = - ^- A + m\ + ^ <pg(£, t) + ^ G 6 (f , £ , t)j <p b (x, t) , 

9b{x, y,t) = - ^— A +m„ + ^ ^(z, i) + ^ , x, t^j 6(x - y) . 
From the variations of the source term K c , we obtain : 

~c ^ <^c ~ c / -* ( \ • $ K c 

w K {x,t) 1 = - — — — , w K {x,t) 2 = t 



5if c (x,t) '"' ' ' ~ 5ir c (x,t) 
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w c 



K 



w c 



5K 

w^'-'msair*- (R13) 

(x*,t) 1 = J*(x,t)+2 J d d x 2 (j^{x,x 2 ,t)^ c {x 2 ,t) + J* u {x,x 2 ,t)n c {x 2 ,t)) , (B.14) 

x, t) 2 = i J u (x , t)+2i J d d x 2 (j* n (£, x 2 , t) p c {x 2 ,t) + J nn (£, x 2 ,t) n c (x 2 ,t)) , 

(B.15) 

l c K (x,y,t) n ^-J^(x,y,t) , (B.16) 

T K (x,y,t) 12 = -2i J^x^t) , (B.17) 

T c K (x,y,t) 22 = J uu (x,y,t) . (B.18) 

The expressions for the matrices t, T, r and R which appear in the dynamical equations 
for the two-time correlation functions are the following : 

tn(x,t) = -g (0 \x,x,t) , t 22 = -1 , (B.19) 

T 1M (x,y,t) = b -^\x,t)5(x - y) , T 2Jk = , (B.20) 

rn 1 i(f,y,t) = -^^(f,*)(J(f-y) , (B.21) 

R 11 , 11 (x,y,t) = ^5(x-y) . (B.22) 
The other matrix elements are equal to zero. 
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